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Abstract 

We describe new conformal field theories based on symplectic fermions that can be extrapo- 
lated between 2 and 4 dimensions. The critical exponents depend continuously on the number of 
components N of the fermions and the dimension D. In the context of anti-ferromagnetism, the 
N = 2 theory is proposed to describe a deconfined quantum critical spin liquid corresponding to a 
transition between a Neel ordered phase and a VBS-like phase. 
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I. INTRODUCTION AND SUMMARY OF RESULTS 



This paper concerns some new renormalization group (RG) fixed points inD = d+l = 3 
dimensional quantum field theory based on symplectic fermions. In Id there exists a vast 
assortment of critical points described by 2D conformal field theory [l]. However in 2d 
the known critical theories are up to now comparatively rare. Important examples are 
the Wilson-Fisher fixed points which are known to describe phase transitions in classical 
statistical mechanics in 3 spacial dimensions as a function of temperature 0, [J. Further 
progress in some important 2c? condensed matter systems, in particular superconductivity 
in the cuprates, has been hindered by the lack of understood critical points in 3D. 

This lengthy introduction will serve to summarize our main results. Since the motivation 
for our model initially came from quantum anti-ferromagnets in 2d, we begin by reviewing the 
aspects of this problem that provide some perspective on our work. For a detailed account 
of this subject with additional references to the original works see Fradkin's book 4|. 

The Heisenberg hamiltonian for a collection of spins S{ on a d dimensional lattice is 

H = JJ2 S t -S 3 (1) 

<i,3> 

For the anti-ferromagetic case ( J > 0), a continuum limit built over a staggered configuration 
close to the Neel state leads to the euclidean action 

S=^Jd D x (d^n-d^n) +S e (2) 

where n 2 = 1 and d^d^ = X^Li ■ The latter constraint makes it an 0(3) non-linear 
sigma-model. In any dimension, the topological term S$ arises directly in the map to 
the continuum when one formulates the functional integral for the spins Si using coher- 
ent states^, 0], and is related to the area swept out by the vector n on the 2-sphere. In Id 
it is given explicitly by 

So= 7T~ I d 2 X EavU ■ {daft x d u n ) (3) 
2tt J 

with 9 = 2ns where s is the spin of S: Sf = s(s + 1). When s is an integer, Sg has no 
effect in the functional integral: the model has no non-trivial infra-red (IR) fixed point and 
the model is gapped. For half-integer spin chains, Sg =7T modifies the IR behavior and the 
model has an infra-red fixed point which necessarily has massless degrees of freedom. This 
is the well-known Haldane conjecture^]. 



The above understanding of Id spin chains relies to a considerable extent on the Bethe- 
ansatz solution 8] , which shows massless excitations. What is more subtle, and was only 
realized much later, is that the low-lying excitations are actually spin 1/2 particles referred 
to as spinons |9[]. These spinons are responsible for destroying the long-range Neel order. The 
n field is composed of two spinons, so this is a Id precedent to the 2d deconfined quantum 
critical points discussed by Senthil et. al.[l(|. In the latter work, numerous arguments 
were given that in 2d there should exist novel critical points that describe for example the 
transition between Neel order and a valence-bond solid (VBS) phase, and this idea strongly 
motivated our work initially. 

In 2c? the term Sg does not appear to have a significant role. One way to anticipate 
this is that unlike in Id where ft and Sg are classically dimensionless (in fact Sg is exactly 
marginal), in 2c? ft has classical dimension 1/2 so that Sg is already RG irrelevant before 
any quantum anomalous corrections. 

Remarkably in 2d a non-trivial IR fixed point appears in the Heisenberg anti-ferromagnet 
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that doesn't rely on the existence of Sg and can be understood in the following way 
12] . The non-linear constraint n 2 = 1 renders the non-linear sigma model perturbatively 
non-renormalizable in 2d. (Unlike in Id, see If a fixed point is understandable by 

Wilsonian RG, this non-renormalizability is potentially a serious problem. However the 
infra-red behavior is captured by the following scalar field theory: 

Swf = J d 3 x Q^n • d^n + A {ft -n) 2 ^ (4) 

where now ft is not constrained to be a unit vector and hence is a linear sigma model 
but with interactions. The above model can be studied in the epsilon expansion around 
D = 4 and the IR fixed point is seen perturbatively. (See for instance fl^ . Ilf]].) This is 
the universality class of the Wilson-Fisher (WF) fixed point, even though it is a quantum 
critical point [16| at zero temperature. We will refer to this model as the O(M) linear sigma 
model and the WF fixed point conformal field theory as in D < 4 dimensions. 

A large part of the literature devoted to the on-going search for other ground states of 
quantum spins represents the ft field in terms of spinon fields z: 

ft = z^az (5) 

where a are the Pauli matrices and z = (zi, z 2 ) = {z{\ is a two component complex bosonic 
spinor. The constraint n 2 = 1 then follows from the constraint z*z = 1. Coupling z to a 



£7(1) gauge field with the covariant derivative 7J M = <9 M — iA^, then by eliminating the 
non-dynamical gauge field using it's equations of motion, one can show that the following 
actions are equivalent: 

J d D x \d^n ■ d^n = J d D x |Z^z| 2 (6) 

(One needs <% • a H = 2S u S jk - SijS H .) 

Since the spinon fields have classical dimension 1/2, in terms of the spinon fields, the 
(n -n) 2 term in the action (jlj) is a dimension 4 operator which is irrelevant in 2d. In analogy 
with the WF fixed point, in order to deal with the non-renormalizability it is natural to 
relax the constraint n 2 = \z^z\ = 1 and to consider (z^z) 2 terms of dimension 2 which are 
RG relevant. However the fixed point is then still in the universality class of the WF fixed 



point. In an effort to perturb the WF fixed point, the authors in [lCj make the gauge field 
dynamical by adding {F^ v ) 2 and thus consider a QED-like theory (in 3-D). It was proposed 
that the model has a fixed point in the universality class of a hedgehog suppressed 0(3) 
sigma model, however because of the expected non-perturbative nature of the fixed point, 
it hasn't been possible to compute any of it's critical exponents. It has also proven difficult 
to see such a second-order transition in simulations of the model [1 71] . 

With the above background we now present the central idea of this paper. The WF 
fixed point in this anti-ferromagnetic context is a quantum critical confined phase since it 
is described in terms of the n degrees of freedom. A deconfined quantum critical point is 
defined then as one describable with the spinon z degrees of freedom. We will postulate 
that for a deconfined critical point the spinon field z should actually be a fermion field, 
henceforth denoted x, and described by the action 

S x = J d D x (d llX %X + A |x f x| 2 ) (7) 

where \ is a two-component complex field, x^X = Tlii=i2x\xi- The non-linear constraint 
X^X = 1 is obviously relaxed. Note that our model contains no gauge field. In terms of real 
fields each component can be written as % = Vi + ^72? X* = Vi ~ ^72 and the free action is 

'S'symplectic = i J d X €ij B^Tjid^Tjj (8) 

where e\2 = —€21 = 1. The 2N real component version has the symmetry 77 — > Urj where U 
is a 27V x 2N dimensional matrix satisfying U T e^U = where = e <8> In- This implies 
the theory has Sp(2N) symmetry, hence is sometimes referred to as a symplectic fermion. 



This kind of theory has a number of important app 



ications in Id, for instance to dense 
19| | . It is known to be a non- unitary 



polymers [18], and to disordered Dirac fermions in 2d 
theory and this potential difficulty will be addressed below. 

Arguments suggesting that symplectic fermions are natural in this context are the fol- 
lowing. First of all, in Id the spinon is neither a boson nor a fermion but a semion, i.e. half 
fermion as far as its statistics is concerned, so there is a precedent for this kind of modified 
statistics. We remind the reader that there is no spin-statistics theorem for 2d relativistic 
theories since spin is not necessarily quantized. Note also that the identity ([6]) remains true 
if z is a fermion. Secondly, suppose the theory is asymptotically free in the ultra-violet. 
Then in this free, conformally invariant limit, one would hope that the description in terms 
of n or z are somehow equivalent, or at least have the same numbers of degrees of freedom. 
One way to count these degrees of freedom is to study the theory at finite temperature 
T = 1/(3 and consider the free energy. For a single species of free massless particle the free 
energy density is 

^S/(S log(lTe "* E) (9) 

— * 

where uj^ = \k\ and the upper/lower sign corresponds to bosons/fermions. In 2d 

^=-c 3 ^T 3 (10) 

where C3 = 1 for a boson and 3/4 for a fermion. (In Id the analog of the above is T = 
— c7rT 2 /6, where c is the Virasoro central charge 2CJ, |2l|.) Therefore one sees that the 3 
bosonic degrees of freedom of an n field has the same C3 as an N = 2 component x field. 
This simple observation is what first pointed us in the direction of symplectic fermions. It 
suggests a kind of bosonization where 3 bosons are equivalent to 4 fermions. In Id, one boson 
is equivalent to 2 fermions, and one can explicitly construct the fermion fields in terms of 
bosons, but we won't need to attempt the analog here. 

Lastly, and most importantly, the symplectic fermion theory has an infra-red stable fixed 
point that is not in the WF universality class, and this is the main subject of this paper. The 
exponents can be computed in the very low order epsilon expansion around D = 4 and they 
are in excellent agreement with the exponents found numerically for the hedgehog-free model 
studied by Motrunich and Vishwanath [22| . The agreement is better than we anticipated. 
We find 

77 = 3/4, i/ = 4/5, ,3 = 7/10 (N = 2,D = 3) (11) 



compared to v = .8 ± 0.1, (3/v = .85 ± 0.05. The shift down to 3/4 from the classical 
value T] — 1 is entirely due to the fermionic nature of the x fields. Since we are discussing 
a zero temperature critical point, temperature exponents are probably not important for 
comparison with experiments, so we also define a magnetic exponent 5 so that (n) ~ B 1 / 5 
where B is an external magnetic field and find 5 = 17/7. It is the fermionic statistics of x 
that is also ultimately responsible for the largeness of these exponents in comparison with 
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fixed point. We thus conjecture that our 2-component model describes the 
jwhich is believed to be a deconfmed critical point of the kind discussed 



the bosonic W 
fixed point in 
by Senthil et. al. 
our model. 

Let us return now to the issue of topological terms. So far we have ignored the gauge 



10]. It is important to point out that there are no "emergent photons" in 



field in the action (Q. Wilczek and Zee have shown how to include a topological term 
that is intrinsic to 2d, which is a Chern-Simons or Hopf term: 

5 cs = J d 3 x e^F^ (12) 

It is well-known that i? modifies the statistics of the z, and this provides an obvious mecha- 
nism for obtaining the fermionic statistics of the x field when *& = ir. It should be emphasized 
that this Chern-Simons term has nothing to do with the S$ term discussed above. Whereas 
the latter arises directly in the map to the continuum, the coefficient of the Chern-Simons 
term is more subtle. The consensus is that for a square lattice d = (j. I24 
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26. 



27|. 



The idea that the fermionic x model should correspond to the action ([6]) with the addition 
of a Chern-Simons term with d = ir helps to resolve the problem that the symplectic fermion 
theory is non-unitary, since the Chern-Simons description is thought to be unitary. In this 
description the x particles are to be viewed as having ir flux attached microscopically. We 
return to the issue of the non-unitarity in section IV and give a different possible resolution 
of it based on a simple projection of the Hilbert space onto pairs of particles. 

In this paper we do not address what kind of microscopic theory can give rise to $ 7^ 
and whether it is related for instance to the spin s of the S, as in 9 = 2irs. But we 
can nevertheless state a 2d version of the Haldane conjecture: When 6 = (mod 2n), the 
quantum critical point is confined and in the universality class of the WF fixed point. On the 
other hand when 9 = it (mod 2ir), the quantum critical point is in the different universality 
class of the fixed point in the N = 2 component fermionic theory described in this paper. 
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The bulk of this paper analyzes the RG for the symplectic model to support the above 
statements. In section II we compute the beta function in D = 4 from the effective potential, 
and display the fixed point in D < 4. In section III we define the critical exponents and relate 
them to anomalous dimensions in the symplectic fermion theory. The relevant correlation 
functions are computed to lowest orders in position space and this determines the exponents 
as a function of N and D. (In terms of Feynman diagrams, which we manage to avoid, 
this involves one and two loop diagrams. For the exponent v one actually only needs 1-loop 
results.) 

In section IV we study our model in a hamiltonian framework in momentum space and 
show how the non-unitarity in manifested. We argue that a simple projection onto pairs 
of particles renders the Hilbert space unitary. We also discuss possible applications to 
superconductivity in the cuprates. 



II. RG BETA FUNCTION 



A. Functional RG 

The 3D fixed point can studied systematically in the epsilon expansion around D = 4^, 



151 ]. The Feynman diagram techniques developed for the WF fixed point, namely dimensional 
regularization, lagrangian counterterms, etc, is easily generalized to our fermionic theory. 
However since we will work to lowest orders only, we need not develop the epsilon expansion 
in much detail since the required quantities can be computed directly in 4D. Working with 
position space correlation functions also helps to clarify the physical content. 

For the beta function it turns out to be convenient to work with the Coleman- Weinberg 
effective potential. This avoids Feynman diagrams and can be especially useful if the po- 
tential has complicated anisotropy, though this will not be investigated here. It also helps 
to keep track of the all-important fermionic minus signs. 

Let x denote a vector of complex fermionic (Grassman) fields x — (Xij • • • j Xn) = {Xi\i 
and consider the euclidean action: 

S x = J d D x {2d^d, X + U( X \ X)) (13) 

where U is the potential. For the purposes of Grassman functional integration, it is con- 
venient to arrange Xi X a 2iV dimensional vector \l/ = (x,X )■ Note that = \I/ T E 1 



/() 1 . 

where in block form Si = , so that is not independent of \I/ as far as functional 

\1 

integration is concerned. 

Consider the functional integral 

Z = I Dm e' s (14) 



in a saddle point approximation. Expanding \I/ around \l/o an d performing the fermionic 
gaussian functional using the formula 



/ 



one obtains Z ~ e 5off where the effective action is 



S cS = S(m )-^TrlogA (16) 



and the Trace is functional. The operator A is 



f-d 2 0\ 1 
A — - -U" (17) 

where £/" is a matrix of second derivatives: 

( d 2 U d 2 U \ 

Assuming [/" is constant, the trace can be computed in momentum space. Dividing by the 
volume one defines the effective potential 

— wiM(:i)-H 

Since we are interested only in the term that determines the beta function in D = 4, we 
expand the log to second order in U" which involves l/k A . Performing the integral over k 
with an ultra-violet cut off \i: 

" d 4 k 1 log/* 
(2tt) 4 £; 4 8tt 2 1 J 



one finds 



where in block form 



V eS = U + ^ tr(S 3 £/"S 3 £/") (21) 



1 
-1 



(22) 
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The renormalization group of the potential then follows from requiring 

V (S (U,n) = V< s (U(a),e a fi) (23) 

We adopt the usual convention in statistical physics and define the beta function as the 
flow toward the infra-red, i.e. to low energies: dll/d£ = —dU /da where e e is a length scale. 
Finally one gets the beta function: 

^ = —^tr(Z 3 U"X 3 U") (24) 
dt 128tt 2 v 3 3 ; v ; 

The flow toward the IR then corresponds to increasing I. 

The RG flow equation (124")) is a functional RG since it determines the flow of the entire 
potential. It is then clear that not all potentials are renormalizable: only if \xi^ 3 JJ"H 3 JJ"') 
is proportional to U is the theory closed under RG flows. If additional terms are generated 
in the trace, they must be included in U until one obtains something renormalizable. By 
"renormalizable" we here mean that the flow just amounts to the flow of some couplings. 



B. Beta function and fixed point 

Let us take the following normalization of the coupling A 

U X = 16tt 2 A U, U= |x f x| 2 (25) 
Denoting <9j = d/dxi, dj = d/dx^ one has 

dfrU = -2xJx] (26) 
d]diU = -2%x f X + 2xJxi 

Evaluating the trace in eq. (j24p one finds that it is proportional to U since U was chosen 
to be isotropic. We also need that the classical dimension of x is (D — 2)/2 which implies 
the dimension of A is 4 — D. Since the leading linear term always has a slope equal to the 
classical dimension of A, the beta function is 

^ = (4- J D)A + (iV-4)A 2 (27) 

The above beta function has a zero at 



Note that A* changes sign at N = 4. It is not necessarily a problem to have a fixed point 
at negative A since the particles are fermionic: the energy is not unbounded from below 
because of the Fermi sea. Near A* one has that d\/d£ ~ (D — 4) (A — A*) which implies the 
fixed point is IR stable regardless of the sign of A*, so long as D < 4. When D > 4 we have 
a short distance fixed point that is not asymptotically free. 



III. SCALING THEORY AND CRITICAL EXPONENTS 

There are two aspects of the scaling theory that differ from the usual WF fixed point for 
classical phase transitions. The first is that n is a composite field in terms of the x's- The 
second is that due to the fermionic nature of x, some of the exponents are in fact negative. 



A. Definition of the exponents for the n field 

Though the spinons x are deconfined, it is still physically meaningful to define exponents 
in terms of the original order parameter n , which is represented by 

ft = x^x (29) 

When N = 2, a are the Pauli matrices; for general N they should be taken as representations 
of the Clifford algebra appropriate to spinor representations of O(N), however we will not 
need these details in this paper. We then define the exponent rj as the one characterizing 
the spin-spin correlation function: 

<™( X M(°)>~^1^ (3°) 

For the other exponents we need a measure of the departure from the critical point; these 
are the parameters that are tuned to the critical point in simulations and experiments: 

S X ^S X + J d D x{m 2 X ] X + B-n) (31) 

Above, m is a mass and B the magnetic field. For classical temperature phase transitions, 
typically m 2 oc (T — T c ). The correlation length exponent u, and magnetization exponents 
(3, 5 are then defined by 

i~ m - u , (ft) ~ mP ~ B 1 ^ (32) 
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Above (n) is the one-point function of the field n (x) and is independent of x by the assumed 
translation invariance. 

To streamline the discussion, let [X] denote the scaling dimension of X in energy units, 
including the non- anomalous classical contribution which depends on D. An action S nec- 
essarily has {S} =0. Using [cPx] = —D, the classical dimensions of couplings and fields are 
determined from [S 1 ] = 0. The above exponents are then functions of [n] and [m]. Since 
[£] = —1, one has v = — [£]/[m] = l/[m]. This, together with eq. fl30|) . implies 

r? = 2[n] + 2- D, u = l/[m] (33) 

One also has 

P = [ft]/[m] = ^(D-2 + T,) (34) 

The magnetic exponent is treated similarly. Treating B as a coupling, then [£>] + [n] = D, 
which implies 

g= D-jn] = D ± 2- 1 

[n] D + rj-2 y 1 

B. Relation to % field exponents 

The fundamental exponents of the symplectic fermion theory are the anomalous dimen- 
sions 7 X , 7 m of x an d m, defined as: 

\x] = (D-2)/2 + lx , [m] = l-j m (36) 
The 7 X exponent determines the two point function of the x fields: 

(X f (x)x(0)) ~ ^ D \ +2lx , (37) 
whereas v is completely determined by 7 m : 

- = 1 - 7m (38) 

The scaling dimension [n] is not a simple function of 7 X to all orders since n must be 
treated as a composite operator. However, since n is quadratic in Xi let us assume that to 
lowest order \n \ =2\x\- With this assumption one has 



r] = D - 2 + 4 7x (39) 
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C. Computation of correlation functions 

To compute 7 x ,7 m we need to consider the following correlation functions in AD. Intro- 
duce an ultra-violet cut-off \x as follows J d 4 x. — > <i 4 x. Suppose the one-point function 
of x}x satisfies 

(x f x) ~ m 2 (l - 2a log n) « (m/i" a ) 2 (40) 
Then the ultra-violet divergence is canceled by letting m — > m(/i) = ra/t Q . This implies 

7m = a (41) 

This is equivalent to the statement \x i x\ = D — 2 + 2^ m . 

The exponent 7 X is determined from the two-point function. Suppose 

(x!(x)x,(0)> ~ ^(1 - 2a'log|x|) « (42) 

Then 

7x = <*' (43) 
We now describe the lowest order contributions to the needed correlation functions. One 
needs the free 2-point functions: 

<x!(x)x,-(0)) = -<Xi(x)x}(0)> = (when A = 0) (44) 

O /I -A. 

Expanding the functional integral perturbatively in A, to first order one has: 

<X f X>A = (x f x)o - 16tt 2 A J d'y (x t x(0)f/(y))) (45) 

Using Wick's theorem, the integrand goes as l/|y| 4 and the ^-integration gives log//: 

(X ] x)x = (l - A(l - N) log//) ( X t x)o (46) 

From dimensional analysis (x*x)a=o oc m 2 . Thus 

= ^ (47) 

The negative sign in 1 — iV in 7 m is entirely statistical in origin. 
To compute 7 X we need to go to second order: 

(x!(x) X ,(0)) A2 = ±(16vr 2 A) 2 / rf 4 y / rf 4 z ( X l(x) X ,(0) C/(y) C/(z)) 



^• il 32^- /(|X|) (48) 
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where J(|x|) is the integral: 

J(|x|) = [d'y [ d'z — * (49) 

J J |y| 2 |x - z| 2 |y - z| b 

It is evaluated in Appendix A where it is shown: 

T(\ \ \ 27T 4 log |x| 

7 (l x l) = r^p (50) 

This in turn implies: 

(1 - N)\* 

7x = 1 ( 51 ) 

D. Values of exponents for arbitrary N, D. 

At the fixed point, one substitutes A = A* into the above expressions. The x field 
exponents are: 

(4 -£>)(! -AT) (4 - Df{l - N) 

7m ~ 2(4 - N) ' 7x ~ A{A-NY 1 ' 

which in turn imply the following n exponents: 

2(4 - N) 



\ - — 1/ / . \ ' LJ -|- 4 

P = 



(4 - N)(D(1 — N) + 2(N + 2)) 

The rj, 5 exponents are given in terms of (3, v in eqs. f l34ll35|) . 

For N = 2, the above expressions give the results quoted in the introduction: v = 
4/5, (3 = 7/10, 7] = 3/4 and 5 = 17/7. 

IV. HAMILTON! AN DESCRIPTION AND UNITARY PROJECTION 

In this section we give a hamiltonian description of our model in momentum space and 
address the non-unitarity issue in a different way than in the Introduction. We also suggest 
how our work may provide some clues toward the understanding of superconductivity in 



the cuprates, which is believed to be a 2 + 1 dimensional problem [3 lj. To do this, one 
must turn to the language of the Hubbard model. In the anti-ferromagnetic phase of the 
Hubbard model, the spin field n = c^ac, where c are the physical electrons. Therefore in 
applying our model to the Hubbard model, the symplectic fermion x is a descendant of 
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the electron, so it seems it can carry electric charge. Consider the zero temperature phase 
diagram of the cuprates as a function of the density of holes. At low density there is an 
anti-ferromagnetic phase. Suppose that the first quantum critical point is a transition from 
a Neel ordered to a VBS-like phase and is well described by our symplectic fermion model 
at N = 2. Compelling evidence for a VBS like phase has recently been seen by Davis' 
group 



32|; and it in fact resembles more a "VBS spin glass". The superconducting phase 
actually originates from this VBS-like phase. Beyond the first deconfined quantum critical 
point at higher density, it is then possible that the 2-component x fields capture the correct 
degrees of freedom for the description of this VBS-like phase. It is important to point out 
that we are imagining that the spinons are still deconfined in the VBS-like phase, in contrast 
to ideas in [h]] . These fermionic spinon quasi-particles acquire a gap away from the critical 
point, which is described by the mass term in eq. fl3TT) . Note that away from the quantum 
critical point, the particles already have a gap m because of the relativistic nature of the 
symplectic fermion. 

Particles with a relativistic kinetic energy are actually not entirely new in 2d physics: 



graphene appears to have massless particles described by the Dirac equation[29|, |30], rather 
than the symplectic fermion theory we considered. In graphene the origin of the massless 
Dirac equation is simply band theory on a hexagonal lattice, so the origin of massless rel- 
ativistic particles is entirely different (and much simpler) than the origin of our symplectic 
fermions. 

The hamiltonian of the symplectic fermion is 

H = y rf 2 x (2dt X %X + 2Vx f • V X + mVx + A (x 1 *) 2 ) (54) 
Expand the field in terms of creation/annihilation operators as follows 

-^(a k e- kx + 6 k e jkx ) (55) 



and similarly for x > where c<j k = Vk 2 + m 2 . Canonical quantization of the x-fields leads to 

{b[, b k >} = -{a k , a k /} = <5 kjk / (56) 

The free hamiltonian is then 

H = J d 2 k u; k ^a k a k + & k & k ) (57) 
14 



Because of the minus sign in eq. (|56|) this is a two-band theory: 



H 6+10) = uj k bl\0) (58) 
H 4|0) = -w k 4|0) 



A two-band structure of this kind has been observed experimentally 32|] and in the Hubbard 



model 



33|]. 



The unconventional minus sign in the anti-commutator of the a, a) is a manifestation of 
the non-unit arity. In particular it implies that the state |k)„ = 4|0) has negative norm: 
(k|k') a = — 5kk'- On the other hand, the two-particle state |ki, k 2 ) a has positive norm. We 
thus propose to resolve the non-unitarity problem by simply projecting the Hilbert space 
onto even numbers of a-particles. It is clear that if they arise from deconfinement of the ft 
field, they will be created in pairs. 

Let us turn to the interactions. In the VBS-like phase the x-particles are charged fermions 
and it's possible that additional phonon interactions, or even the x 4 interactions that led 
to the critical theory, could lead to a pairing interaction that causes them to condense into 
Cooper pairs just as in the usual BCS theory. Recent numerical work on the Hubbard model 



suggests that the Hubbard interactions themselves can provide a pairing mechanism [33]. 
The (x^X) 2 interaction is very short ranged since it corresponds to a 5-function potential 
in position space. Because of the relativistic nature of the fields, the interaction gives rise 
to a variety of pairing interactions. Let us examine pairing interactions within each band 
that resemble BCS pairing. If all momenta have roughly the same magnitude |k|, then the 
interaction gives the term (up to factors of n): 

H in t = -A ^2 ( a k,i a -kj a -k,i«k,i + (a -> b) J + .... (59) 

k; i,j=U 

The overall minus sign of the interaction is due to the fermionic statistics. To compare with 
the BCS theory, the interaction contains terms such as 

H in t = ( a kT a -ki a -ktaki) + •••• (60) 

k 

Because of the overall minus sign this is an attractive pairing interaction as in BCS. One 
difference is that some pairing interactions have the spins flipped in com par ison with BCS. 



This discussion has been expanded and a few more results derived in [341 ]. 
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V. CONCLUSIONS 



We have described new 3D fixed points based on symplectic fermions and have applied 
the N = 2 case to deconfined quantum criticality. In this interpretation, the symplectic 
fermions are the deconfined spinons. The main evidence for our model is the agreement 
with the exponents found in [22j for the hedgehog-free 0(3) sigma model. The fermionic 
nature of the x-field is what is responsible for the negative contributions to exponents like 
rj that bring it below the classical value r\ — 1. 

Our model is non-unitary, however we have addressed this in two ways, one based on a 
unitary Chern-Simons description, the other based on a projection of the Hilbert space onto 
even numbers of particles. 

After this work was completed, new simulations by Sandvik report evidence for a decon- 
fined quantum critical point in a Heisenberg model with four-spin interactions. There it was 



22j and with our prediction of v = 4/5. The 



found that v = 0.78 ± 0.03 consistent with 
r) exponent on the other hand, rj = 0.26 ± 0.03 is quite different from both our result and 



the one in 22j. This could simply mean that the critical point in the four-spin interaction 
model is in a different universality class. It could also mean there are significant corrections 
to our exponents at higher order or due to the compositeness of n ; recall the fixed point at 
lowest order occurs at the relatively large value A* = 1/2. Higher order computations are 
currently in progress and we hope to report our results in the near future. 

Our exponents are well defined for any iV < 4, including negative N. Based on the 
comparison of exponents, our model for N a negative integer was proposed to describe 



O(M) models such as the Ising model in 36] with the identification M = —N 
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VII. APPENDIX A 



In this appendix we do the integral I(|x|) in eq. (l4"9j) . Shifting z = z' + y, and using the 
identity 



1 rl 



AB J (tA + (1 - t)B) 2 

one can do the integral over y: 

■L 



dt 7Ta — T~7i 77"m? (61) 



/ ^Vr^ - = 7r 2 (log(|w| 2 /L 2 )-3) (62) 

Jo |y||y-w| 2 



where w = x — z' and L is an IR cut-off. Introducing an ultra-violet cut-off /i, J <i 4 x 



J^li (i 4 x the integral over z' can be performed giving 



7T 



1 



J(|x|) = — (2 log(|x|/L) + 1) - 6ttV (63) 
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